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Optical Rogue Waves in Vortex Turbulence 


Christopher J. Gibson, Alison M. Yao, and Gian-Luca Oppo 

SUPA and Department of Physies, University of Strathelyde, Glasgow G4 ONG, Seotland, U.K. 

We present a spatio-temporal mechanism for producing 2D optical rogue waves in the presence of 
a turbulent state with creation, interaction and annihilation of optical vortices. Spatially periodic 
structures with bound phase lose stability to phase unbound turbulent states in complex Ginzburg- 
Landau and Swift-Hohenberg models with external driving. When the pumping is high and the 
external driving is low, synchronized oscillations are unstable and lead to spatio-temporal vortex- 
mediated turbulence with high excursions in amplitude. Nonlinear amplification leads to rogue waves 
close to turbulent optical vortices, where the amplitude tends to zero, and to probability density 
functions (PDFs) with long tails typical of extreme optical events. 


Rogue waves (RWs) in high seas were once a thing of 
legend: massive walls of water tens of metres high, ca¬ 
pable of destroying large ships, appearing from nowhere 
then vanishing leaving no trace of their existence fl. Be- 
cause of the very short lifetime of such extreme events, 
quantitative studies and simulations of the mechanisms 
behind their creation in oceanography have grown only 
recently ii- Although the origin of these waves is still 
under debate, RWs have been realized in several optical 
systems Q from optical fibres (HI, to optical cavities 
and photonic crystals 

Here we investigate a spatio-temporal mechanism in¬ 
volving vortices in a 2D turbulent state that is capable 
of generating RWs, building upon previous work con¬ 
cerning a singly resonant optical parametric oscillator 
system with a low amplitude detuned seeding field j§|. 
For generality and application to a variety of nonlin¬ 
ear dynamical systems, we consider a Forced Complex 
Ginzburg-Landau (FCGL) equation 0 and a Forced 
Complex Swift-Hohenberg (FCSH) equation, both with 
external driving. We focus on the loss of synchronisa¬ 
tion of the Adler locked states obtained at large driving 
amplitudes. When decreasing the external forcing, oscil¬ 
lations at the Adler frequency become spatio-temporally 
unstable leading first to a phase and then to an amplitude 
instability that forces, locally and randomly, the forma¬ 
tion of pairs of oppositely charged vortices. Since the 
total power in the transverse direction remains almost 
constant throughout, the nonlinearity pushes the inten¬ 
sity to high spikes close to interacting vortices, resulting 
in the rare formation of RWs. The RWs described here 
are outside thermodynamic equilibrium, do not survive 
in the purely temporal (single mode) case and are due 
to a deterministic, nonlinear and vortex-mediated turbu¬ 
lence far removed from a purely stochastic superposition 
of optical waves. 

To demonstrate the generality of optical RWs in vortex 
turbulence we employ a variety of mathematical models: 

dtE = Ein -{l-iuj)E-^ iV^E + P f{\Ef) E 

-V{uPeV^fE, ( 1 ) 

where E is the complex field, Ejn is the (real) amplitude 
of the external forcing, uj is the frequency difference be¬ 


tween the unperturbed field and the external driver, 
is the transverse Laplacian, P is the laser pump, f{\E\‘^) 
is 1 — \E\‘^/3 for the laser [Tl| and sinc^{\E\) for the 
optical parametric oscillator (9||. Time is normalised to 
the photon decay rate in the optical cavity and space to 
^LX/Att where L is the cavity length and A the wave¬ 
length. Finally, F is zero for the FCGL model and one 
for the FCSH case, where e is a small parameter (here 
fixed at 0.01) due to the fast dynamics of atomic vari¬ 
ables in lasers [l2|. The FCGL and FCSH models can be 
applied in many othe r sy stems, e.g. chemical oscillations 
(l3| . granular media [l3| and hydrodynamics 0- 

The cases of relevance are obtained when the detun¬ 
ing UJ is different from zero. In this case the frequency 
locked states that one observes at large driving ampli¬ 
tudes become unstable upon decreasing the driving Ejn. 
For fixed values of uo and P, the homogeneous stationary 
states of Eq. m have a typical P-shaped dependence 
on ^IN as displayed, for example, in Fig. [1] where the 
stability of these solutions to perturbations of zero wave- 
vectors is shown. The uppermost lines in the P-shaped 
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FIG. 1. (Golor online) Stationary intensity of plane waves 
and their stability (stable = solid green, unstable with real 
eigenvalues = dashed blue, unstable with complex eigenval¬ 
ues = dash-dotted red) for the FGGL model with P = 4 and 
(jj — 0.53. The black dotted lines represent the minima and 
maxima of stationary hexagonal patterns, the vertical lines 
where the optical turbulent state starts (shaded area). The 
blue circle is the stationary intensity of the laser with no in¬ 
jection. The FGSH model displays very similar results. 


curves of Fig. [T] correspond to the homogeneous locked 
states where the external driving is large enough to over¬ 
come the frequency difference with the injected device. 
When increasing P/at, a saddle-node bifurcation heralds 
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the onset of the frequency and phase locked homogeneous 
states. When, instead, decreasing the parameter the 
homogeneous solution loses stability to spatially periodic 
patterns with a critical wave-vector given by kc = 

In Fig. [T]the maximum and minimum intensities of the 
hexagonal patterns obtained numerically when re¬ 
ducing the external driver are displayed via a black dot¬ 
ted line. Although the phase of the pattern is periodically 
modulated in space, the stationary character of these pat¬ 
tern solutions demonstrates that they are locked to the 
frequency of the injection. A typical hexagonal structure 
in the case of the FCGL equation with finite size input 
beams is shown in Fig. Note that all the results pre¬ 
sented in this paper remain valid in the limit of transverse 
periodic boundary conditions. As the amplitude Ej^ of 
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FIG. 2. (Color online) Transverse intensity for (a) hexago¬ 
nal Turing pattern and (b)-(c) optical vortex-mediated tur¬ 
bulence. Corresponding transverse phase distribution of two 
oppositely charged vortices (d). Simulations of the FCGL 
equation with parameters, P = 6,^; = 0.77, Ein = 1.00 (a) 
and Ein = 0.95 (b)-(d), see [3l. The beam radius is IOtt. 
(c) and (d) correspond to the area of the white square in (b). 

the external drive is further reduced (see shaded area in 
Fig. m, spatially periodic patterns become unstable and 
a regime of unlocked dynamics sets in i[i3- Figure [3] 
shows the temporal evolution of an unstable hexagonal 
pattern in an Argand {Im{E) versus Re{E)) diagram. 
The hexagonal pattern (see Fig. is a phase bound so¬ 
lution that progressively loses stability along a circle in 
the Argand diagram corresponding to a phase instability 
(Fig. [SJd). This phase instability then grows into an am¬ 
plitude instability (Fig. [3t) that leads to the formation 
and annihilation of pairs of oppositely charged vortices 
(see Fig. [Jt and d) and a regime of spatio-temporal ir¬ 
regularity similar to the vortex-mediated turbulence de¬ 
scribed in in the CGL in the absence of forcing (see 
Fig. [31i). The helical waves propagating around the de¬ 
fects act as the driving force behind the turbulent state. 



A typical instantaneous intensity distribution of this tur¬ 
bulent state is presented in Fig. [2]3 & c. The interacting 
vortices correspond to the localized regions of zero ampli¬ 
tude (shown in black). Note that the turbulent dynamics 
of vortices is deterministically driven by the spatially cou¬ 
pled nonlinearity and not by the superposition of random 
waves typical of optical speckles [18|. Indeed, in the case 
of speckle, the field distribution in the Argand plane has 
a Gaussian shape as opposed to the almost circular one 
shown in Fig. [3li. There are also noticeable intrinsic dif¬ 
ferences in the field correlations 0 and in the PDFs of 
the intensity (see below). 
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FIG. 3. (Color online) Field distributions in the Argand plane 
of the unstable hexagonal pattern at t=0 (a), during phase 
instability (t=183) (b), in a regime of amplitude instability 
(t=236) (c) and in a turbulent state (t=472). The green cir¬ 
cle in (d) is the Adler limit cycle. Simulations of the FCGL 
equation with the same parameters as Fig. [T]and Ein = 0.60. 

To better understand the nature of the turbulent state 
in forced models CD we consider dynamical solutions cor¬ 
responding to unlocked oscillations and their robustness 
to spatially dependent perturbations in the FCGL model 
(F = 0) where analytical predictions are feasible. In the 
absence of spatial coupling an approximate limit cycle 
trajectory for the field can be found by period averaging 
methods, E = Aq [cos(0(t)) + isin(0(t))], [HI, where 
Aq = 3(P — 1)/P, 0(t) is well approximated by its period 
average Qt tt with Q = \JuR — and ujl = Ejn/Aq. 
When Q is real, i.e. in the absence of locked states, the 
trajectory is the phase-drift solution of the Adler equa¬ 
tion [HI, dtcj) = uo — uol . Such solution is clearly 

phase unbound and is superimposed onto the Argand di¬ 
agram in the turbulent regime of Fig. [3li to show that its 
underlying dynamics is ruled by the unlocked state. The 
accuracy of the approximate solution has been checked 
for a wide range of Ejn values in the FCGL model. The 
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excursions in intensity do not exceed 10% while those in 
frequency are well within 1%. We have then proceeded 
to study the stability of the spatially synchronized oscil¬ 
lation by including spatial coupling in the FCGL. The 
stability eigenvalues of the spatially synchronized limit 
cycle are given by 

A± = -{P - 1) ± ^(P - 1)2 -{to- fc2)2 (2) 


where k is the spatial wave-vector. At the critical wave- 
vector for pattern formation kc = ftie stability A+ is 
marginal but perturbations due to the approximate na¬ 
ture of E induce a slow instability of the synchronous 
oscillation. The eigenvector associated with A+ is along 
the limit cycle, again demonstrating a phase instability. 
As mentioned earlier, this phase instability grows into 
an amplitude instability and then into vortex-mediated 
turbulence as demonstrated numerically in Fig. IH start¬ 
ing from low amplitude noise. A homogeneous zero state 
with added noise quickly evolves towards the unstable 
limit cycle (from 0 to 4 in Fig. [4^). The limit cycle 
dynamics first synchronises the spatial oscillations (see 
the narrow line at t = 90 in Fig. [ 4 ) 3 ) and then moves 
towards the vortex turbulence state via phase (Fig. IDd) 
and amplitude (see Fig. [3]:) instabilities. We outline 
that the mechanism of spontaneous vortex creation in 
the FCGL and FCSH is not trivial. In contrast with 
the CGL, stationary vortex solutions are not possible in 
driven systems like o as all locked states have bound 
phases around that of the injection. However, at low 
driving amplitudes, moving vortices and vortex-mediated 
turbulence in ([T]) are possible due to the Adler unlocked 
dynamics of the limit cycle trajectory. It is known [22 1 
that the adiabatic elimination of the polarization variable 
introduces an all wave-vector instability of the spatially 
homogeneous state below the point where the linear sta¬ 
bility of the lower branch of the ^'-shaped homogeneous 
state predicts complex conjugate eigenvalues (see Fig. 
[ 1 ]). This feature, in principle, may have serious conse¬ 
quences in the turbulent regimes. A second important 
consequence of our analysis, however, is that Eq. CD 
for F = 0 displays a very fast dynamics that takes the 
system towards the limit cycle where large wave-vector 
instabilities are promptly eliminated (see Fig. [4])). The 
large wave-vector instability of the lower branch of the 
homogeneous stationary states is not present in the case 
of FCSH when F = 1. 

This mechanism for vortex turbulence is essential for 
the generation of RWs in externally driven systems de¬ 
scribed by the spatio-temporal dynamics of Eq. CD- 
These systems are outside thermodynamic equilibrium, 
do not display relaxation oscillations and present a deli¬ 
cate balance between the energy input and the losses (the 
pump P, injection P/w, nonlinearity /(|Pp) and output 
mirror). During the turbulent evolution, the total power 
V remains almost constant at values close to those of the 
laser with no injection (see Eig. [5^). By considering the 
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FIG. 4. (Color online) Field distributions in the Argand plane 
for the FCGL equation with the same parameters as Fig. [1] 
and Ein = 0.70. (a): t=0 (black), t=0.40 (red), t=0.75 

(blue) and t=1.4 (green), (b): t=90 (black), t=256 (red), 
t=280 (blue) and t=314 (green). 


energy density and the energy flux of the FCGL equation 
the time evolution of the power is given by 


^ = 2/ 

dt J 


Eij^Re^E) 




dxdy 


(3) 

where (x, y) is the transverse plane. For the approximate 
limit cycle solution the power V is conserved at the value 
of ttwqAq where wq is the beam width of the input laser. 
In the turbulent state, however, maintaining an almost 
constant power in the presence of moving vortices of zero 
intensity implies the simultaneous appearance of large 
amplitude spikes. If the vortex density is large, multi¬ 
vortex collisions can occur with the production of large, 
short-lived spikes in the intensity (see Fig. [Sb). Short¬ 
lived large intensity spikes are rare but possible events, 
fitting the characteristics of RWs. The particular shape 
and symmetry of these spikes is crucially determined by 
the number and position of the surrounding vortices, a 
feature that is unique to this particular mechanism of 
RW formation. RWs in single transverse mode class-B 
lasers with injected signals have been observed in Q but 
due to relaxation oscillations and not to 2D vortex tur¬ 
bulence. In fact, without spatial coupling due to diffrac¬ 
tion, no RWs can be observed in systems described by 
Eq. ([1]). To characterize our spatio-temporal RWs, we 
use a commonly accepted definition of statistically rare 
events SHS.i-e. if the intensity of the field at a spatial 
point over a long period of time is greater than the mean 
wave height plus eight standard deviations then the wave 
can be classified as an extreme event or RW, similarly to 
the significant peak intensity method Q. 

In Eig. [6] we show PDEs for different regimes of vortex 
turbulence. The wave statistics is well-approximated by 
a Gaussian fit when the pump intensity is low (blue solid 
line in Eig. E]). At higher pump intensities the long-tailed 
PDEs show mass generation of RWs. Here the statistics 
changes drastically and is very well approximated by a 
Weibull distribution [3, [l3| (red dashed line in Eig. E]). 
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FIG. 5. (Color online) (a) Time evolution of the power for 
three values of the pump (P = 2, 4 and 5) and for the FCGL 
(solid lines) and FCSH (dashed lines) models during vortex 
turbulence. The horizontal lines are the predictions from the 
approximate Eq. ([3]). (b) RW spike in the transverse inten¬ 
sity for the FCGL with P — — 0.77, P/at = 1.00. The 

peak intensity is 42.22, the average is 4.93 and the standard 
deviation is 1.57 (all in normalised units). 


Note that non-Gaussian PDFs cannot be replicated by 
superpositions of random waves. We also note that the 
RWs in vortex turbulence demonstrated in Fig. [6] are 
different from those due to vorticity in models of inviscid 
fluids (H. The CGL and CSH have been shown to be 
equivalent to the flow of a compressible and viscous fluid 
with density p = | Pp and velocity v=V0 where (j) is the 
phase of the field [26|. In the case of our forced systems, 
Vxv remains extremely close to zero in the locations 
where RWs are observed. We conclude that our RWs are 
due to the interaction of free vortices in the absence of 
vorticity. 



FIG. 6. (Color online) Intensity PDFs for the FCGL model: 
P — 2^uj — 0.3, Ein — 0.24 (purple triangles) with Gaussian 
fit (blue line) and P = 8,Cc; = 2.4, Ein — 3.40 (green stars) 
and the sinc^ model: P — — 1, Ejn — 1.48 (black circles) 

with Weibull fit (red dashed line) [1^. The vertical black line 
is the threshold for defining waves as an extreme event. Inset: 
Parameter region for RWs (red area) in the FCGL. The upper 
curve is the pattern to vortex-mediated turbulence transition, 
the lower curve the turbulence to target pattern transition. 
No RWs are observed below the red dashed line. 

Finally, we show in the Fig. [6] inset the wide parameter 
region where we observe RWs induced by vortex turbu¬ 
lence in systems with external driving of the FCGL kind. 
Very similar results have been obtained for the FCSH 


and different nonlinearities such as those of the FCSH 
and the singly resonant OPO, thus demonstrating the 
universality of the phenomenon. 

In conclusion, we demonstrate a mechanism for pro¬ 
ducing RWs in the transverse area of externally driven 
nonlinear optical devices via vortex turbulence. Given 
the universality of our model, this mechanism should be 
observable in a large variety of systems. Models of lasers 
with injected signal, where the invariance of the Adler 
limit cycle is well known I 111 . l24l | , can be easily extended 
to semiconductor media [2^ and to class B lasers, thus 
including the largest majority of solid state lasers. Out¬ 
side optics, vortex-mediated turbulence without driving 
has been observed in nematic liquid crystals (131 , chemi¬ 
cal reactions jl^l and fluid dynamics [29[. In the unlocked 
regime of these systems with driving, vortex turbulence 
can excite RWs and lead to the formation of highly inho¬ 
mogeneous fields with non-Gaussian statistics. 

The prototype model used to describe RWs is the 
Non-Linear Schrodinger (NLS) equation dSi. The 
FCGL and FCSH models studied here are active, non¬ 
conservative systems outside thermodynamic equilibrium 
where many of the methods developed for the NLS can¬ 
not be applied. In the NLS equation, as well as in the 
CGL and CSH equations, stationary vortex solutions are 
possible although mainly unstable. In the presence of 
forcing, vortices can only exist in dynamical states. It 
is exactly in these situations that we have demonstrated 
RWs close to regions of interaction of turbulent vortices. 
Because of universality, suitably perturbed NLS models 
may also display these features. 

Unlike RWs in the longitudinal direction Q, the aspect 
ratios required for transverse RWs induced by 2D vortex 
turbulence are extremely small (typical input beams have 
diameters less than 1mm) and the statistics require times 
of the order of hundreds of ps. The small aspect ratio, the 
full 2D character and the quick dynamics represent the 
major advantages of transverse optical devices in study¬ 
ing the generation and control of RWs with applications, 
by universality, in hydrodynamics and oceanography. 
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